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Abstract 

We prove the Bochner-Weitzenbock formula for the (nonhnear) Laplacian on gen- 
eral Finsler manifolds and derive Li-Yau type gradient estimates as well as parabolic 
Harnack inequalities. Moreover, we deduce Bakry-Emery gradient estimates. All 
these estimates depend on lower bounds for the weighted flag Ricci tensor. 

1 Introduction 

Geometry and analysis on singular spaces is an important topic of current research. Be- 
sides Alexandrov spaces, Finsler manifolds constitute one of the most relevant classes of 
explicit examples of metric measure spaces. Finsler spaces quite often occur naturally via 
homogenization as scaling limits of discrete or Riemannian structures. 

A Finsler manifold is a smooth manifold M equipped with a norm - or, more generally, 
a Minkowski norm - F{x, ■) on each tangent space T^M. The particular case of a Hilbert 
norm leads to the important subclasses of Riemannian manifolds. In our previous paper 
[OSlj . we introduced and studied in detail the heat flow on a Finsler manifold. It can 
equivalently be defined 

• either as gradient flow in L'^{M,'m) for the energy ^ jj^ F{^uY dm\ 

• or as gradient flow in the L^-Wasserstein space {V2{M), W2) for the relative entropy 
Jj^j. u log u dm. 

For this nonlinear evolution semigroup, we proved C^'"-regularity and L^-contraction 
estimates as well as integrated upper Gaussian estimates a la Davies ( jOSll Sections 3, 
4]). 

Various flne properties of the heat flow are intimately linked with the weighted flag 
Ricci tensor Kicn (with N G [dimM, 00]) introduced in |Ohl] . In |Ohl] . the flrst author 
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proved that the curvature- dimension condition CD{K,N) in the sense of Lott-Sturm- 
Villani for a weighted Finsler space (M, F, m) is equivalent to a lower bound K for RIctv- 
In particular, the relative entropy - regarded as a function on the L^-Wasserstein space - 
is i^'-convex if and only if RiCoo is bounded from below by K. In [OSH Theorem 5.6], we 
proved that lower bounds for the weighted flag Ricci curvature imply point-wise compar- 
ison results a la Cheeger-Yau. Surprisingly enough, however, in |0S2] we observed that 
exponential contraction/expansion bounds in Wasserstein distance 

W2{Ptf^,Ptiy)<e'''W2{fi,iy) (1.1) 
hold true if and only if we are in a Riemannian setting. 

The main goal of this paper now is to derive the Bochner- Weitzenbock formula (The- 
orems 13. 3[ 13. 6p 

A^"(^^^^) -/^(At.)(V«) = RiCoo(V«) + ||V^t.||L(v.)- (1-2) 

Here A^" denotes the linearization of the Laplacian A on M in direction Vm. The 
formula fll.2p a priori holds true only on the set {x G M | V'u(x) ^ 0}. It extends to an 
identity in distributional sense on entire M. 

Formulation and proof of (11. 2p - as well as applications of it - require lot of care due 
to the lack of regularity. For instance, solutions to the heat equation Au = du/dt on M 
will be only if M is Riemannian ( |OSlj ). As a first striking application we deduce the 
Bakry-Emery gradient estimate (Theorem 14.11) 

F{VPtu{x)Y < e~'^'PZ-{F{Vuf){x). (1.3) 

Here Pt denotes the (nonlinear) semigroup generated by the Laplacian A, whereas -Pq^" 
denotes the (linear, symmetric) Markov transition operator on L^(M, m) with time- 
dependent generator L^, < s < t, obtained by linearization of A in direction V-u(s, ■). 
An immediate consequence of (11.31) is the growth bound for Lipschitz constants 

Lip(PiM) < e"^*Lip(M). (1.4) 

Note that in the Riemannian setting (or more abstract linear frameworks), according to 
the Kantorovich-Rubinstein duality, the bounds (11. ip and (ll.4p are equivalent to each 
other (with C = -K). 

Finally, we prove the Li- Yau gradient estimate (Theorem 14. 4p for positive solutions to 
the heat equation on Finsler spaces which in the case Riciv > simply reads as 

9 N 

F(V(logn)(t,x)) -a<(logn)(t,x) < -, 
and deduce as a corollary the famous Li- Yau type Harnack inequality (Theorem 14. 5p 

ft\^''^ f dix yf\ 
u{s^x)<u{t^y).[-^ exp^^^J 

for any < s < t and x,y E M. 

For (11.31) and the subsequent estimates, M is required to be compact. It is an open 
question whether this compactness assumption can be replaced by completeness. 



2 



2 Finsler geometry 



We briefly review the fundamentals of Finsler geometry, for which we refer to |BCS] and 
[Shj . as well as some results from |Ohl] and |0S1] . 



2.1 Finsler manifolds 

Let M be a connected, n-dimensional C°°-manifold without boundary. Given a local 
coordinate (a;*)"^^ on an open set U C M, we will always use the coordinate {x'^,v^)^^i of 
TU such that 

" d 

V = y^v' — 

i=l 



e TrM, X eU. 



Definition 2.1 (Finsler structures) We say that a nonnegative function F : TM — ¥ 
[0, oo) is a C°° -Finsler structure of M if the following three conditions hold: 

(1) (Regularity) F is on TM \ 0, where stands for the zero section. 

(2) (Positive 1-homogeneity) It holds F{cv) = cF[v) for all v G TM and c > 0. 

(3) (Strong convexity) The n x n matrix 

is positive-definite for all v G T^M \ 0. 

We call such a pair (M, F) a C°° -Finsler manifold. 

In other words, F provides a Minkowski norm on each tangent space which varies 
smoothly in the horizontal direction. For x,y ^ M, we define the distance from a; to ?/ in 
a natural way by 



d{x,y) :-- 



inf / F(T]{t))dt, 
^ Jo 



where the infimum is taken over all C^-curves r] : [0, 1] — > M such that ?7(0) = x and 
ri{l) = y. Note that d{y, x) ^ d{x, y) may happen since F is only positively homogeneous. 
A C°°-curve 77 on M is called a geodesic if it is locally minimizing and has a constant speed 
(i.e., F{fj) is constant). See (12. 6 p below for the precise geodesic equation. For v G T^M, 
if there is a geodesic rj : [0, 1] — > M with fj{0) = v, then we define the exponential map 
by exp^{v) := ''^(l). We say that {M,F) is forward complete if the exponential map is 
defined on whole TM. Then the Hopf-Rinow theorem ensures that any pair of points is 
connected by a minimal geodesic (cf. |BCSl Theorem 6.6.1]). 

For each v G T^M \ 0, the positive-definite matrix (fi'ij('i^))"j=i in (12.11) induces the 
Riemannian structure of T^-M via 

9^[J2^^d^.lJ2b^d^ J - J:9^Av)a.b,. (2.2) 
^ 1=1 j=i ' «j=i 
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This is regarded as the best Riemannian approximation of F\t^m in the direction v. In 
fact, the unit sphere of is tangent to that of F\t^m at v /F{y) up to the second order. 
In particular, we have gy{v,v) = F{vY. The Cartan tensor 

is a quantity appearing only in the Finsler context. Indeed, Aijk vanishes everywhere on 
TM \ if and only if F comes from a Riemannian metric. We will repeatedly use the 
next useful fact on homogeneous functions. 

Theorem 2.2 (cf. |BCS[ Theorem 1.2.1]) Suppose that a differentiable function if : \ 
— y M satisfies H{cv) = d'Hiv) for some r G M and all c> and f G M" \ {in other 
words, H is positively r- homogeneous) . Then we have 

= rHiv) 

1=1 

for allv G M"\0. 

Observe that gij is positively 0-homogeneous on each T^M, and hence 

n n n 

Ai,j,{vy = A,,k{vy = Y = (2-3) 

i=l j=l k=l 

for all V G TM \ and i, j, k = 1,2, . . . , n. Define the formal Christoffel symbol 

1=1 ^ ■' 

for V G TM \ 0, where {g'^^) stands for the inverse matrix of (gij)- We also introduce the 
geodesic spray coefficients and the nonlinear connection 

" 1 flC^ 

G\v) := Y l}kivyv', N]{v) := -^{v) 

j,k=i ^ 

for V G TM \ 0, and ^^(O) = A/'j(0) := by convention. Note that G' is positively 
2- homogeneous, so that Theorem 12.21 implies Yl]j=i-^ji'^)'^'' ~ G'*(f). 

By using the nonlinear connections Nj, the coefficients of the Chern connection is 
given by 

■■= l]k - E -f^^3imK + Afc^AT- - A,,^iV-) (2.4) 

l,m=l 

on TM \ 0. Namely, the corresponding covariant derivative of a vector field X = 
J2i=i X^{d/dx'^) by f G T^M with reference vector w G T^M \ is defined by 

d:x{x) := E ^^5^(^) + E^}.(-K^'(^)LTL• (2-5) 

i,j=l ^ k=l ^ ^ 
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Then the geodesic equation is written as, with the help of (12. 3p . 

_d_ 

dx 



Dm = E {^^(^) + G'm)}^ _ = 0. (2.6) 

i=l 



2.2 Nonlinear Laplacian and the associated heat flow 

Let us denote by C* : T*M — y TM the Legendre transform. Precisely, C* is sending 
a e T*M to the unique element v e T^M such that a{v) = F*{aY and F{v) = F*{a), 
where F* stands for the dual norm of F. Note that C*\t*m becomes a linear operator 
only when F\t^m is an inner product. For a differentiable function u : M — > M, the 
gradient vector of m at x is defined as the Legendre transform of the derivative of u, 
Vu{x) := C*{Du{x)) e T^M. If Du{x) = 0, then clearly Vu{x) = 0. If Du{x) ^ 0, then 
we can write in coordinates 

We must be careful when Du{x) = 0, because gij{Vu{x)) is not defined as well as the 
Legendre transform C* being only continuous at the zero section. For later convenience 
we set in general 

Mv:={xeM\ V{x) ^ 0} (2.7) 

for a vector field V on M, and M„ := Mvu. For a differentiable vector field V on M and 
X G My, we define VV"(a;) G T*M ® T^M by using the covariant derivative (12.51) as 

VV{v) := D^V^x) G T.M, v G T^M. (2.8) 

We also set V^'u(x) := V(V'u)(x) for a twice differentiable function u : M — t- M and 
X G M„. 

Lemma 2.3 'V^u{x) is symmetric in the sense that g'^u{^'^u{v),w) = g-^u{^'^u{u)),v) 
for all f , w G TxM. 

Proof. This can be checked by hand. Indeed, choosing a coordinate {{d/dx^)\x}^=i or- 
thonormal with respect to g-vu{x), we see by calculation 

j=l k=l ■' 



dx^ 



□ 



From now on, we fix an arbitrary positive C°°-measure m on M as our base measure. 
Define the divergence of a differentiable vector field on M with respect to m by 
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where we decompose m in coordinates as dm = e* dx^dx^ ■ ■ ■ dx"^. This can be rewritten 
(and extended to weakly differentiable vector fields) in the weak form as 



/ (f)diVmV dm = — / D(f)(V) dm 
J M J M 



for all (f) G C^(M). Then we define the distributional Laplacian of m G Hl^^{M) by 
Au := diVmCVu) in the weak sense that 



■pAudm := — / D(j){'Vu) dm 
M J M 

for G C^(M). We remark that Hl^^{M) (and Lf^^^M) etc.) is defined solely in terms of 
the differentiable structure of M. 

Given a vector field V such that 7^ on M„, we can define the gradient vector and 
the Laplacian on the weighted Riemannian manifold {M,gv,m) by 

on M \ M„, 

where the latter should be read in the sense of distributions. Note that, in the definition 
of A^u, we used the divergence with respect to m rather than the volume form of gy- It 
is not difficult to see V^"u = Vu and A^"u = A.u ( |OSH Lemma 2.4]). We also observe 
for later use that, given u, fi and /2, 

W2(V^"/i) = ^7v«(V^"/i, V^"/2) = I^/i(V^"/2). (2.9) 

In [OSlj . we have studied the associated nonlinear heat equation du/dt = Am. The 
nonlinearity is inherited from the Legendre transform. Given an open set Q C M, define 
the Dirichlet energy of m G Hl^^{Q) by 

S^(u) := I I FiVuYdm. 



We will suppress f2 if = M, namely S = £m- Set 

H\n) := {u G Hl^{n)fM'^{yi,m)\En{u) < 00}, 
and let HliVl) be the closure of C'^{VL) with respect to the (Minkowski) norm : = 

Definition 2.4 (Global and local solutions) (1) For T > 0, we say that a function 
u on [0,T] X M is a global solution to the heat equation if u G L^([0, T], iJQ(M)) fl 
H^{[Q,TIH--^{M)) and if 



I 4>-K~ dm = — f D(j){^Ut) dm 
holds for all t G [0,T] and G C^{M), where we set Ut := u{t, 
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(2) Given an open interval J C M and an open set Q C M, we say that a function 
M on J X is a local solution to the heat equation on I x Q if u G L'fod^ ^ 



Global solutions can be constructed as the gradient flow of the energy functional 8 in 
L'^{M,m). We summarize the existence and regularity properties established in [OSlj in 
the next theorem. 

Theorem 2.5 f jOSTl Sections 3,4]) 

(i) For each Uq G Hq{M) and T > 0, there exists a global solution u to the heat equation 
on [0, T] X M , and the distributional Laplacian is absolutely continuous with 
respect to m for all t G (0,T). // M is compact, then such a global solution is 
uniquely determined by its initial datum uq. 

(ii) Given an open set Q C M, the continuous version of any local solution u to the heat 
equation on enjoys the Hf^^-regularity in x as well as the 0^'°" -regularity in both t 
andx. Furthermore, the distributional time derivative dtu lies in HI^^{VL) nC(fi). 

We remark that the mild smoothness assumption |USH (4.4)] clearly holds true for 
our C°°-smooth F and m. 

2.3 Ricci curvature 

The Ricci curvature (as the trace of the flag curvature) on a Finsler manifold is defined 
by using the Chern connection (and is in fact independent of the choice of connection). 
Instead of giving a precise definition in coordinates, here we explain a useful interpretation 
due to Shen |Sht §6.2]. Given a unit vector v G T^M (i.e., F{y) = 1), we extend it to 
a C°°-vector field on a neighborhood of x in such a way that every integral curve of 
V is geodesic, and consider the Riemannian structure gy induced from fl2.2p . Then the 
Ricci curvature Ric(w) of v with respect to F coincides with the Ricci curvature of v with 
respect to gv (in particular, it is independent of the choice of V). 

Inspired by the above interpretation of the Ricci curvature, the weighted Ricci curva- 
ture for (M, F, m) is introduced in |Ohl] as follows. 

Definition 2.6 (Weighted Ricci curvature) Given a unit vector v G T^M, let rj : 
{—s,e) — )■ M be the geodesic such that 77(0) = v. We decompose m as m = e~* vol^ 
along rj, where vol.^ is the volume form of (7^. Define 



F{Vu) G Li 



loc 



(/ X Q) and if 




holds for all (p G C^{I x Q). 




(2) Riciv(i;) 



Ric{v) + (^ o r/)"(0) 



(^077)^0)' 
N -n 



for G {n, 00) 
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(3) RiCoo(t^) := Ric{v) + o r]y\0). 

For c > and G [n, oo], we define RicN{cv) := c^RicAr(f). 

It is established in |Ohlt Theorem 1.2] that, for i^' e M, the bound RicAr(t;) > KF{vY 
is equivalent to Lott, Villani and the second author's curvature- dimension condition 
QD{K,N). This extends the corresponding result on (weighted) Riemannian manifolds 
(due to [vRSJ . [StlJ . [St2] . |St3] . |LV2] . |LV3] ). and has many analytic and geometric 
applications (see |Ohl] . |0S1] or a survey |Qh2] ). 

Remark 2.7 (a) In contrast to A^^-u = Am, RicAr(VM) may not coincide with the 
weighted Ricci curvature Ric^"(V'u) of the weighted Riemannian manifold {M^ g-^u-,^)- 
Indeed, on a Minkowski space (M", || ■ ||), it is easy to see that (]R",5fv«) is not fiat for 
some function u. 

(b) For a Riemannian manifold (M, volg) endowed with the Riemannian volume 
measure, clearly we have = and hence RIcat = Ric for all G [n, oo]. It is also 
known that, for Finsler manifolds of Berwald type (i.e., Ff^ is constant on each T^M), the 
Busemann-Hausdorff measure satisfies (\I' o r^)' = (in other words, Shen's S-curvature 
vanishes, see |Sht §7.3]). In general, however, there may not exist any measure with 
vanishing S-curvature (see |0h3j for such an example). This is a reason why we begin 
with an arbitrary measure m. 

For later convenience, we introduce the following notations. 

Definition 2.8 (Reverse Finsler structure) Define the reverse Finsler structure F 
of F by F (v) := F{—v). We say that F is reversible if F = F. We will put an arrow ^ 
on those quantities associated with for example, d{x,y) = d{y,x), = — V(— m) 
and Sic7v(f) = RicAr(— f). 

We say that (M, F) is backward complete if (M, "W) is forward complete. Compact 
Finsler manifolds are both forward and backward complete. 



3 Bochner-Weitzenbock formula 

In this section, we prove the Bochner-Weitzenbock formula. For the sake of simplicity, 
we first derive a general formula for vector fields V, and then apply it to gradient vector 
fields of functions. We start with a point- wise calculation on My (recall 02. 7p ) followed 
by an integrated formula. 

3.1 Point-wise calculation 

We follow the argumentation in [Vl] Chapter 14] in the Riemannian situation, with the 
help of [OhT] and [UST] . 

As in the previous section, let (M, F) be an n-dimensional C°°-Finsler manifold and 
m be a positive C°°-measure on M. Fix a C°°- vector field V on M and x G My. For 
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t G (0, 6] with sufficiently small 5 > 0, we introduce the map Tt and the vector field Vt on 
a neighborhood of x by 

Ttiy) := exp^[tV{y)], V,{rt{y)) := -^{y). (3.1) 

As the curve a{t) := Tt{y) is geodesic, taking the covariant differentiation of t i — 
Vt{a{t)) yields the pressureless Euler equation 



^ + D^V, = 0, (3.2) 



where we abbreviated as 

1=1 1=1 

Now we put rj{t) := Tt{x), take an orthonormal basis {cj}"^]^ of {TxM.gv) such that 
&n = v{^)/ ^{v{^))j consider the vector field along r] 

Ei{t) := D{Tt)M) e T^{t)M, z = 1, 2, . . . , n. 

Then each Ei is a Jacobi field along rj. (Note that En{t) ^ r]{t) / F{'q{t)) in general, since 
Tt{r]{s)) ^ rj{s + 1) for s,t > 0.) We set i?-(t) := D^^Eiit) for simplicity. Define the n x n 
matrix-valued function B{t) by E[{t) = X]j=i (This B corresponds to U in 

|Vi[ Chapter 14].) Along the discussion in |Ohl[ Lemma 7.2], we obtain the following 
Riccati type equation. 

Lemma 3.1 (Riccati type equation) For rj and B as above, we have 
d{trB) 



dt 



-(t) +tr(5(t)2) +Ric(r7(t)) = 0, t > 0. 



Proof. We give only an outline of the proof, see |Ohl] for details. Consider the matrix- 
valued function A{t) = {gr,{Ei{t), Ej{t))) and observe A' = BA + AB'^ by definition (5^ 
is the transpose of B). Since each Ei is a Jacobi field, we have [BA — AB^)' = 0, so that 

BA - AB^ = B{0) - B{Of, A' = 2BA - B{0) + B{Of. 

We also deduce from the Jacobi equation for Ei that 

A" = -2 Ric^ +2BAB^ = -2 Ric^ +2B^A - 2B (fi(0) - B{Of) , 

where (Ric^),^ := grji^K^ {Ei, i))r], Ej) and is the curvature tensor. Comparing this with 
A" = 2B'A + 2BA', we find B' = - Ric^ A'^ - B^. Taking the trace with respect to Qr, 
completes the proof. □ 

Lemma 3.2 (i) We have B(t) = 'Wt{ri{t)) in the sense that, for each i = 1,2, . . . ,n, 

vVt{E,{t)) = Y.UhAt)EAt)- 
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(ii) It holds that tr(i?(t)) = divm Vtivit)) + where m = e * vol^ along r] such 

that vol^ denotes the Riemannian volume measure of Qr). 

Proof, (i) By the definition dH]), VVt{Ei{t)) = L'^*(^)\4(?7(t)). Setting 



Ei{t) 



d_ 



0=0 



fc=i 



we have 



/c=l j,k,l=l 

where rj^t) = Yl^=iV''(t){9/dx^)\rj(t)- Exchanging the order of differentiations (by S and 
t) in the first term, we deduce from (13.11) that 



^(*) = |K-(7;(exp.&.)); 



<5=0 



Therefore we obtain E'-{t) = D^^^^Vt{ri{t)) = 'Wt{Ei{t)) and complete the proof. 

(ii) Choose a coordinate (x*)"^]^ around ri{t) such that {{d / dx^)\ri(t)}i'=i is an orthonor- 
mal basis of (T^^t^M, gvt). We will suppress evaluations at rj{t). Recall first that 



d_ 



k=l ^ 



d 

dx^ 



Thus we have 



A \dvl 



tr{B{t)) = tr{VV,) = E ^ + J2^mV^ 

i=l ^ j=l 

We also find, by ([MD and (1^ . 



Next we observe from dm = e *-\/det^i4 (ix^dx^ ■ ■ ■ fix" on rj that 



Note that 



a(log v^det gvt) _ 1 / dgvt 



dx^ 



2" 



-1 



dx"- 



■9v, 



j=i fc=i 
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This yields, together with (13. ip . 



yi djlog ^/det gy,) ^ 1 



^(v;)v7 + %f(v<)'f'(«) 



l,j=l 



dv 
F{Vt) 



G\Vt) 



We used the geodesic equation f l2.6p of t] in the second hne. This completes the proof. □ 
We deduce from Lemmas I3.2( ii) that 
d 



dt 



div^\4(r?(t)) + tr (5(0)2) +Ricoo(?7(0)) = 0. 



Thanks to (13. 2p . we have 

d 
It 



t=o+ 



dV 

dwmVt{r^{t)) = Z}(div„ V)(77(0)) +div„ ( ^ 



i=0+ 



Combining these, we obtain 

diY^{DlV) - D{dw^ V){y) = RiCoo(V^) + tr (^(O)^) 



(3.3) 



at X. 



Now, we take u G C°°(M) and x G M„, and apply (13. 3p toV = Vn. Then the symme- 
try of V^ti (Lemma 12. 3p allows us to simplify (13.30 in two respects. First, Lemma [3.2( i) 
immediately yields tr(_B(0)2) = || V2m(x)||^^(^^), where || ■ ||hs'(Vm) stands for the Hilbert- 
Schmidt norm with respect to gvu- Second, for each i = 1, 2, . . . , n, 

d ( 5'v«(Vm, Vm) 



This implies 



3\ 9 _„^jF(Vu 



Plugging these into (13.31) . we obtain (13. 4p below on M„. 

Theorem 3.3 (Point-wise Bochner-Weitzenbock formula) Given u G C°°(M), we 
have 

(3.4) 
(3.5) 



, Vu 



{^^^) - D{Au){Vu) = Ric^(V«) + \\V\\\HSi^u) 



as well as 



A 



2 



- D{Au){Vu) > Ricjv(VM) + 



[An) 

N 



for N G [n, oo], point-wise on Mu 
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Proof. We have seen that (13 ■4p holds on M^. The second assertion is clear if iV = oo. 
For N e (n, oo), we observe from Lemma E^l that, as {ci}^^^ is orthonormal so that -B(O) 
is symmetric. 



(tr5(0))^ 



> 



(Am + L)^(Vm)) 



n 



Plugging a = Am and b = D\E'(Vm) into 



{a + bf 



a 
N 



we obtain 



n I\ N -n 

(Am + L)^(Vm))2 



+ 



N{N - n) 



n 



> 



(Au) 



5(0)- 



tr5(0) 



n 



HS 



a b 
N ^ N-n 



n N N-n 

Together with (13 ■4p . this yields the desired estimate (I3.5p . The remaining case of N = n 
is derived as the limit. □ 



Remark 3.4 Our Bochner-Weitzenbock formula (13. 4p can not be simply derived from the 
formula of the weighted Riemannian manifold {M, g-^uy'm)- This is because RiCoo(V'u) 
and II V^mII^^,^^^-) are different from those for (yfvu (unless all integral curves of Vm are 
geodesic), while D{Au){'Vu) and A'^"(F(V'u)^/2) are same for F and g-^u- Recall Re- 
mark [2]Tl^a) and note that V^'u(t') ^ (V'^'^)^u{v) in general (for v ^ Vm). 

Let us close the subsection with additional comments on the situation of \^ = Vu 
with u G C^{M). A little knowledge of optimal transport theory is necessary for which 
we refer to ^i] and |Ohl] . A potential function (pt of the vector field Vt (i.e., Vy^t = Vt) 
is given by 

d{x,yf 



inf 



2t 



+ u{x) >. 



(3.6) 



This is the Hopf-Lax formula providing a (viscosity) solution to the Hamilton- J acohi 
equation 

that corresponds to (13.21) . Indeed, we can apply the proof of [LVH Theorem 2.5] verbatim 
to see (13. 7p . although our distance is nonsymmetric. In view of optimal transport theory, 
(13. 6p is rewritten as 



inf 



d{x,yf 



{-tu{x))\ = {-tur{y), 



where {—tuY stands for the c-transform of the function —tu for the cost c{x,y) 
d{x,yY/2. Thus, if we put rjit) = exp^lfVu^x)], then it holds that 

f,{t) = -^(-<^,)(r^(t)) = Vipt{rj{t)). 
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Hence, in terms of Wasserstein geometry, the vector field Vy^t is the tangent vector at time 
t of the geodesic (exp[tV'u]jj/x)ig[o,5] starting from any absolutely continuous probability 
measure fi (whereas 6 depends on u), where exp[)f:VM]jj/i denotes the push- forward of /x by 
the map exp[tV?i]. 

3.2 Integrated formula 

We need to be careful in extending the point-wise estimate fl3.5l) to a global one (in the 
weak sense), because some quantities are undefined on M \ M„. We prove an auxiliary 
lemma to overcome this difficulty. 

Lemma 3.5 Let L be a locally uniformly elliptic, linear second order differential operator 
in divergence form on a Riemannian manifold. Let£^ denote the associated Dirichlet form 
with domain Hq{M) and with intrinsic gradient denoted by V^, i.e., 

g^{h,h) = - [ h-Lhdm= [ \V^h\'^dm, h e H^{M). 

J M JaI 

Then we have the following. 

(i) For each h e Hl{M), it holds V^h = a.e. on h^^{Qi). 

(ii) If he Hl{M) n L°°(M), then V^{h^/2) = hV^h = a.e. on /^-^(O). 

(iii) The assertions in (i) and (ii) also hold true if h merely lies locally in the respective 
spaces. 

Proof, (i) is a particular case of the general fact that, for local regular Dirichlet forms 
with square field operator F, 

E^{h,h)= [ T{h,h)dm= I T{h,h)dm. 

(In the monograph |BHj . this is called the energy image density.) 

(ii) follows from the chain rule. 

(iii) Given an arbitrary open relatively compact set U C M, choose a cut-off function 
e Hl{M) n L°°(M) such that = 1 on t/. Then, for each h e i^ioc(^), it holds 
hep e H^iM) and we have V^(/i0) = a.e. on (hcpy^O) nU = 0-^(0) n U. Thus we 
can reduce the claim to (i), and it is similar for (ii). □ 

Theorem 3.6 (Integrated Bochner-Weitzenbock formula) Given u E H^^^(M) fl 
C^{M) with Am e Hl^{M), we have 

= [ <p{D{Au){Vu) + RicUVu) + \\V'ufHSiVu)}dm 
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for all G H^iM) f] L°^{M), and 

for N G [n, oo] and all nonnegative functions (p £ H^{M) fl L°°{M). 



dm (3.8) 



Proof. As the proofs are common, we consider only fl3.8p with N < oc. 

(I) Let us first treat the case of u G C°°(M). If G Hl{Mu), then follows from 
flS.Sp via integration by parts. 

For an arbitrary bounded nonneg ative function G H^iM) n L°°(M), set 

0fc := min{0,A;2F(VM)2}, A; G N. 
Then 0^ G H^{Mu) and limfc_j.oo 0fc(a^) = ^(2;) for all x G M^. Hence we have 

- j^D(l)^{v^-(^l^^^y^ dm > l^^<Pk\^D{Au){Vu) + RicMiVu) + ^-^^ dm. 

(3.9) 

In the limit as k goes to infinity, the right hand side converges to 

/ (f)I.D{Au){Vu) + RicN{Vu) +^^^\ dm. 

The first two terms of the integrand obviously vanish on the set M \ M^. In order to 
see that also the third term vanishes, let us fix a local coordinate (x*)"^]^ and apply 
Lemma [3.51 to the function h := du/dx^ G Hl^^{M). It implies that dh/dx^ = a.e. on 
h~^{0) for all j = 1, 2, . . . , n. In other words, d'^u/dx^dx^ = for all i, j = 1,2, . . . ,n on 
HiLii^'"/^^* ~ 0}- Hence, in particular. Am = a.e. on M \ M„. 

To see the passage to the limit on the left hand side of (13. 9p . we observe from 
Lemma [33] that V^''(F(Vm)2) = 2F(Vm) V^"(F(Vm)) vanishes a.e. on M \ M„. Thus, 
putting 

Qfc := {x G M„ I 0(x) > Pf(Vm(x))'} = {a; G MJ 0(x) 7^ 0fc(x)}, 

we find 



D(</. - (/.fc) (v^"(F(Vn)2)) rfm 



< 





+ e [ 




Jn 



The first term of the right hand side tends to zero since Qj. decreases to a null set as k 
goes to infinity. For the second term, note that 

D(F(Vn)2) (y^^{F{Vuy)^ = 4F{Vuy ■ D{F{Vu)) (v^"(F(Vm))) 
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on M„. Therefore, by the choice of we obtain 

J D{F{Vuf){v^''{F{VuY)^dm<Aj (j) ■ D{F{Vu)){v^'' [F{Vu))^ dm 

(/c oo). 

Hence, as k goes to infinity, the left hand side of (13 .Qp converges to 

This proves the claim for u G C°°(M). 

(II) Now let us consider the general case of m G Hf^^{M)r^C^{M) with Au G Hl^^{M). 
Choose Uk G C°°{M) with Uk ^ u locally in the if^-norm as A; — cxd. Using integration 
by parts, we rewrite (13. 8p for the function Uk in the form 



M 



2 



>- / div„(0VMfc)AMfc(im + / 0<J RicAr(Vnfe) + -^^^ }> rfm. 



As A; goes to infinity, each of the terms appearing in this inequality converges towards the 
repective term with u in the place of Uk- This finally proves the claim for general u. □ 



4 Applications 

In linear semigroup theory, the estimates as in Theorems 13. 3 j are also called the 
inequality or the Bakry-Emery curvature- dimension condition (after |BEj ). Applying it 
to solutions to the heat equation, we obtain Bakry-Emery and Li-Yau type estimates (for 
N = oo and N < oo, respectively). We remark that Theorem 12.51 ensures that solutions 
to the heat equation enjoy enough regularity for applying Theorem 13.61 Throughout 
the section, we assume the compactness of M. We will consider only global solutions 
u : [0,T] X M — y M to the heat equation (more precisely, their continuous versions on 
(0,T] X M) for simplicity. 

Given a global solution u : [0,T]xM — > M to the heat equation, let us fix a measurable 
one-parameter family of non- vanishing vector fields {V^}rG[o,T] on M such that Vr = 
on Mu^. For < s < t < T, denote by P^^"^ the (linear, symmetric) Markov transition 
operator on L'^{M,m) associated with the time-dependent generator A^. Precisely, for 
each h G L'^{M, m), ht = P^t^'h for t G [s, T] is the weak solution to dht/dt = A^'ht with 
hs = h. 

Theorem 4.1 (Bakry-Emery gradient estimate) Let {M,F,m) be compact and sat- 
isfy RiCoo > K for some G M, and let u : [0,T] x M — > M. be a global solution to the 
heat equation. Then we have 

F{Vut{x)y < e-2^(*-^)p,3"(F(VM,)2)(x) (4.1) 

for alio < s <t <T and x e M. 
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Proof. For fixed t G (0,T], u as above and for an arbitrary nonnegative function h G 
C(M), we set 

P^-h■^^^^dm, 0<s<t. 



His) :-- 



^2Ks 



Then we deduce from the definition of -P,"^" that 



H'(s) 



-,2Ks 



D 



M 



Note that, thanks to Lemma [3. 5 [ the integrals in the right hand side are well-defined. By 
(12. 9p . the first term in the right hand side coincides with 



dm. 



The other terms are calculated as, with the help of Theorem 12.2 

' du 



^2Ks 



Pj-h{ D 



'M 

^2Ks 



ds 



dm 



' j P3''^|d(Am,)(Vu,) + KFiyusf^ dm. 
To be precise, we used the positive 0-homogeneity of gij to see, on M^^ 







d ' 


-1 




ds 



1 -A i.,— ^dusdus 



n 



dx^ dsdx^ ' 2 ^ dv'' ''^ 

i,j=l i,j,k=l 



D 



D 



dus 
ds 



ds dx^ dx^ 
d{Vus)k 



ds 



(4.2) 



Therefore we have H'{s) < by the weak formulation of the Bochner-Weitzenbock formula 
(13. 8p (for (f) = PYt'h and = oo) together with the hypothesis RiCoo > K. 
Thus H is non-increasing and, in particular. 



^2Kt 



h ^ dm < e^^' 



M 



PsTh ■ 



dm 



M 



^2Ks 



h-P. 



Vm 

s,t 



M 



FiVusf 



dm. 



Since this holds true for any nonnegative h, we obtain the claimed inequality (14. ip for a.e. 
X G M. The Holder continuity in x of both sides of (14. ip finally allows to deduce it for 
every x G M. □ 
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We stress that it is used in (14. ip the mixture of the nonhnear operator Ug i — > Ut and 
the hnear one P'^i^- As an immediate corollary to Theorem I4.H we obtain a growth bound 
for Lipschitz constants. For a continuous function u G C(M), define 

Lip(m) := sup — ^ — . 

x^y(zM d[x,y) 

Corollary 4.2 Assume that (M, F, m) is compact and satisfies RiCoo > K for some K G 
M, and let m : [0, T] x M — Y R &e a global solution to the heat equation. Then we have 

\\F{Vut)\\Loo <e~'''\\F{Vu^)\\L^ 

and, ifuQ G C{M), 

Up{ut) < e"-^^*Lip(Mo) 

for allt G [0,T]. 

Remark 4.3 In the Riemannian case, each of the properties in Corollary 14.21 character- 
izes manifolds with RiCoo > K ( [vRSl Theorem 1.3]), and RiCoo > -ft' is also equivalent to 
the i^-contraction property of heat flow in the L^-Wasserstein space ( [vRSl Corollary 1.4], 
see fll.ip ). The equivalence between the gradient estimate and the contraction property is 
known to hold for more general linear semigroups by Kuwada's duality [Ku] . Comparing 
Theorem 14.11 and Corollary 14.21 with the non-contraction property of heat flow estab- 
lished in [US2j . we observe that a completely different phenomenon occurs for nonlinear 
semigroups. 

Next we consider Li-Yau type estimates for < oo (cf. |LY] and |Da] ) . 

Theorem 4.4 (Li-Yau gradient estimate) Let us assume that {M,F,m) is compact 
and satisfies Kicn > K for some N G [n, oo) and i^' G M, put K' := mm{K, 0} and take 
a positive global solution u : [0,T] x M — > R to the heat equation. Then we have 



F{Vilogu){t,x)Y -9dtilogu){t,x) < NO^ 
on (0,T] X M for any 6 > 1. 



1 K' 



2t A{e - 1) 



Proof. We follow the argumentation in |Dat Lemmas 5.3.2, 5.3.3] and focus on the modifi- 
cations required because of nonlinear ity and the lack of higher order regularity. Through- 
out the proof, we fix a measurable one-parameter family of non-vanishing vector fields 
{Vr}r£io,T] as dcscribcd in the paragraph preceding to Theorem 14.11 We divide the proof 
into five steps. Note that it suffices to show the claim for t = T. 

(I) First consider the function / := logM which is in space and C^'" in both space 
and time. As u solves the heat equation, / satisfies 

Af, + FiVfty = dJt (4.3) 

for every t in the weak sense that 

{-D<j){Vft) + <j)F{Vft)^}dm= [ <f)dtftdm 

M JM 



17 



for each G H^{M). Note also that g-^f^ = g-^ut ^-^^ ^ft ^ H^{M) for each 

t. 

(II) Next we verify that the function a{t,x) := tdtf{t,x) (g H^{M) by Theorem 12.51) 
satisfies 

A^a-9t(7+y + 2Da(V/) = (4.4) 
in the sense of distributions on (0,T) x M. Indeed, for each G ifgllO)^) ^ 

j D(t){V^(T) + (9^0 ■ a + ^ + 20D(t( V/) I cimcit 

{ - D(t0)(V^(9t/)) + dt{t(t>)dtf + 2i0D(aj)(V/)} dmdt 

{D(5i(t0)) (V/) + (A/ + F{Vff) + t09,(F(V/)2) } rfmrft 



JM 

T 





= 0, 

where we used fl4.3p and fl4.2p in the second equality. To be precise, a calculation similar 
to 04.21) ensures 

dt{D{t<p){Vf)) = D{dtm){Vf) + D{t<p){V''{dJ)). 

(Ill) Now let us consider the function a{t,x) := t{F(V/(t, x))^ — 6dtf{t,x)}. It lies 
in (M) for each t and is Holder continuous in both space and time. Using the previous 
identity 04.41) for a as well as our version of the Bochner-Weitzenbock formula, we shall 
deduce that 

A^a + 2Da{Vf) - dta > (3 (4.5) 
again in the distributional sense on (0, T) x M , where /3 denotes the continuous function 

■■= -J + 2t|l(F(V/)2 - djY + KF(V/)^|. 

Indeed, for each nonnegative G H^{{0,T) x M), (|0]), and (jM]) show 

^ / { - ^<^(V^a) + 20Da(V/) + 9*0 ■ a + ^| dmcit 

tZ}0(v^"(F(V/)2)) +2t0D(F(V/)2)(V/) 

JM ^ ^ ^ 

- 0{F(V/)2 + 2tD(9i/)(V/)} + 0F(V/)2} dmrft 

/ { - tD0(v^"(F(V/)2)) - 2t0D(A/)(V/)| M. 
Thanks to 03.81) and the hypothesis RIcat > this is estimated from below by 



JM 



2t I / 0<| KFiVf Y + 1^ rfmrft 
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and (14. 3p completes the proof of (14.51) . 

(IV) Let {s,x) denote the maximizer of a on [0, T] x M (note that M is compact). 
Without restriction, a{s, x) > and thus s > (otherwise, the assertion of the theorem is 
obvious). We shall claim that this implies f3{s,x) < 0. Assume the contrary, (3{s,x) > 0. 
It would imply /3 > on a neighborhood of {s,x). Hence, according to (14.51) on such 
a neighborhood, the function a would be a strict sub-solution to the linear parabolic 
operator 

div^(V^a) + 2Da{Vf) - dta. 

Therefore, a{s,x) would be strictly less than the supremum of a on the boundary of any 
small parabohc cylinder [s — 5, s] x Bs{x), where Bs{x) := {y G M\d{x^y) < 5}. In 
particular, a could not be maximal at (s,x). 

(V) The conclusion of the previous step is that /3(s, x) < 0. In other words. 



a > 



at the maximum point {s,x) of a. Following the reasoning of |Dat Lemma 5.3.3], we 
conclude 

, , N9^ [ K's 
a(s,x) < — 

and thus, since K' < 0, 



1 



X , , NOW K'T 
Oi\T, y) < a(s, x) < 



2 V 2(^-l)_ 

for all y G M . This completes the proof. □ 
We proceed along the line of |Da[ Theorem 5.3.5] and obtain the Harnack inequality. 

Theorem 4.5 (Li-Yau Harnack inequality) Assume that {M,F,m) is compact and 
satisfies Ricjv > K for some N G [n,oo) and K & W, put K' := minjfT, 0} and take a 
nonnegative global solution u : [0,T] x M — > M to the heat equation. Then we have 

, ^ . ^ /t\''^/' fedix.yf OK'Nit-s) 
u[s,x) < u{t,y) ■ I - exp 



sj "^XAit-s) 4(^-1) 

for any 6>l,0<s<t<T and x, ?/ G M. 

Proof. Replacing u hj u + e if necessary, we may assume without restriction that u is 
positive. Let rji^r) = exp^(rw) for r G [s,t] be the minimal geodesic from x = rj{s) to 
y = r](t) with suitable v G T^M. Then obviously F(?7(r)) = d{x,y)/(t — s) for all r. We 
also put / := logM, 

^ ed{x,yf 9K'N , , w . ^^ 
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Then we have, by Theorem H? 

dMr)=Dfir]) + dtf+^-^ + Q 

It 



Hence we obtain 



which proves the claim. □ 
In both Theorems 14.41 and 14.51 we can choose 6' = 1 if = 0. Thus we have the following. 

Corollary 4.6 Let {M,m,F) be compact with Ricjv > for some N G [n,oo), and take 
a nonnegative global solution m : [0, T] x M — )■ M to the heat equation. Then we have 

( \^ (+ \ (^\^'^ fd{x,y)' 
u[s,x) < u[t,y) • I - exp 



Sj \4:{t-S) 



for any < s < t < T and x,y & M. 
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